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Abstract 

A didactical exposition of the classical problem of the trajectory determination of a 
body, subject to the gravity in a resistant medium, is proposed. Our revisitation 
is aimed at showing a derivation of the problem solution which should be as simple 
as possible from a technical point of view, in order to be grasped even by first-year 
undergraduates. A central role in our analysis is played by the so-called "chain rule" 
for derivatives, which is systematically used to remove the temporal variable from 
Newton's law in order to derive the differential equation of the Cartesian representation 
of the trajectory, with a considerable reduction of the overall mathematical complexity. 
In particular, for a resistant medium exerting a force quadratic with respect to the 
velocity our approach leads to the differential equation of the trajectory, which allows 
its Taylor series expansion to be derived in an easy way. A numerical comparison 
of the polynomial approximants obtained by truncating such series with the solution 
recently proposed through a homotopy analysis is also presented. 



1 Introduction 

The study of the motion of a body under the action of the gravity is one of 
the first applications of the Newton's laws met by a student in a typical first- 
year course of general physics. Since the topic is not particularly exciting to 
present, one is tempted to find unconventional ways to expose the subject in 
order to render it more appealing to the audience. The determination of the 
body trajectory, and of the range in particular, is a good candidate for this; if 
the body is in vacuum the laws of motion assume, in the time domain, simple 
analytical expressions that make the trajectory determination a nearly trivial 
task. On the other hand, when the body moves in air, where the resistance 
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acting on it is customarily modeled as a force anti-parallel to its velocity and 
whose modulus is proportional to the first (for "slow" objects) or to the second 
(for "fast" objects) power of the body speed, the mathematical complexity grows 
so rapidly that even nowadays the study of the motion in a resistant medium 
remains an active research field, as witnessed by the recent literature [TJ [3J SI El 

EH El El 0. 

The aim of the present work is to give a didactical presentation of the tra- 
jectory determination problem as simple as possible from a technical point 
of view, in order to be grasped even by first-year undergraduates. We are 
aware of several nice didactical approaches to this problem appeared in the 
past [HI H21 H31 [IH HH HS1 HH HI] ; nevertheless, we believe that it could be 
possible to further reduce the mathematical level of the presentation to avoid 
resorting to concepts, like for instance the hodograph 19, 20], which will become 
available to the students typically during the second year. Our presentation is 
organized as a sequence of steps of growing complexity: starting from the ideal, 
free- friction resistance linear in the velocity is considered, up to the most 

difficult (and still open) problem of the motion of a body in the presence of a 
quadratic resistant force. For all three cases we propose the same approach 
based on the use of Cartesian coordinates, which naturally occur in describing 
the shape of the body trajectory. Another aspect we are willing to emphasize 
is that the solution of several problems in mechanics often does not require the 
knowledge of the complete solution of Newton's equations in the time domain, 
but rather the overall mathematical complexity of an apparent nontrivial prob- 
lem can be considerably reduced by resorting to some "tricks" with which the 
student should become familiar to as soon as possible!]] Among them the use of 
the so-called "chain rule" for derivatives appears to be a technical instrument of 
invaluable help to allow nontrivial problems to be addressed also at an elemen- 
tary level of treatment To help the student to familiarize with this technique, 
the classical parabolic trajectory in vacuum is re-derived without touching the 
temporal laws of motion, but rather by removing the temporal variable from the 
Newton's law, written in the "natural" Cartesian reference frame, through the 
use of the chain rule. This is done to introduce in the most transparent way to 
the student the approach subsequently employed when a linear resistance force 
is added to the dynamical model and which allows the Cartesian equation of 
the trajectory to be exactly retrieved only through elementary quadratures. In 
the final part of the paper the same methodology is applied to the most realistic 
case of a quadratic resistance force, for which the determination of the analyti- 
cal expression of the trajectory equation still remains an open problem. In this 
case, the proposed approach leads to the differential equation for the trajectory 
in Cartesian coordinates, which, within the small slope approximation, can be 
integrated in an elementary way to retrieve the classical solution provided by 
Lamb [TTJ[33]. Moreover, the differential equation turns out to be particularly 

1 Richard Feynman was used to call the collection of such tricks as the "box of tools" of 
any physics student 12 1 1 . 

2 A celebrated example of use of the chain rule in mechanics has been provided by Arnold 
Sommerfeld in deriving the first Kepler's law without resorting to the energy conservation 1221 . 



2 Trajectory of a body in vacuum 



3 



suitable to be solved by a power series expansion, whose single terms can, in 
principle, be evaluated analytically up to arbitrarily high orders. Finally, the 
same differential equation is used to provide an elementary derivation of the 
hodograph of the motion, thus establishing a sort of changeover toward more 
advanced approaches. 

2 Trajectory of a body in vacuum 

We consider the classical parabolic motion of a point mass m under the action 
of the sole gravity force, so that Newton's second law reduces to 

a = g, (l) 
where a denotes the point acceleration and g the gravity acceleration. In Fig. [1] 




Fig. 1: Geometry of the problem. 

the Cartesian reference frame Oxy used to describe the motion of the point mass 
is depicted, with Vq denoting the initial velocity (with x- and y- components 
vq_ x and wn,y, respectively). With respect to this frame Eq. (p} splits into the 
pair of scalar differential equations 

0, 

(2) 

~9, 

where v x and v y denote the x- and y-component of the velocity, respectively, 
and the dot denotes derivation with respect to the time t. The trajectory of 
the motion is naturally described, in Cartesian coordinates, by the function 




3 Trajectory in the presence of a linear resistance force 



4 



y = y(x) and our aim is to extract the related differential equation without 
solving Eq. ([2]) in the time domain. In doing so, we have to formally remove 
from the equations of motion the temporal variable t. This can be done in a 
very simple way, by using the chain rule for derivates. 

To this end we start from the relation r = v written in Oxy, i.e., 

x = v x , 

(3) 

y = Vy, 

and, after dividing side by side the two equations, we obtain 

which coincides with the slope of the trajectory r\ — tan 8, with 9 being the 
angle between v and the x-axis, with the superscript denoting derivative with 
respect to the spatial variable. To find the differential equation for y(x) it is 
now sufficient to derive both sides of Eq. Q with respect to x and to take the 
chain rule, i.e., 

ax v x at 

into account. On substituting from Eq. §5§ into Eq. ((4]) we have 

d 2 y 1 d ( v y \ v y v x - v x v y ^ 



dx 2 v x dt \v, 
which, on taking Eq. ([2]) into account, leads to 

d2 V = _9_ = 9_ , 7 s 

dx 2 v 2 v 2 x ' 

and where, in the last passage, use has been made of the fact that the x- 
component of the velocity is constant. Equation (|7|) can be solved via elementary 
quadratures; in particular, on imposing the "initial" (i.e., at x — 0) conditions 

2/(0) = 0, 

(8) 

y'(0) = tan O , 



the well known equation for the parabolic projectile trajectory is obtained: 

3 Trajectory in the presence of a linear resistance force 



y(x) = x tan 6»o - — |— x 2 . (9) 

0,x 



Consider now the motion of the point mass in the presence of a resistant force 
proportional to the point velocity, i.e., of the form — bv, with b denoting a 
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suitable parameter. In this case Newton's law can be cast in the following form: 

a = g v, (10) 

T 

where the parameter r — m/b rules the temporal scale of the resulting dynamics. 
This, in particular, suggests the use of "natural" and dimensionless normalized 
quantities: r will become the temporal unity, the productor (the so-called limit 
speed) the speed unity, and gr 2 the displacement unityo With such units the 
subsequent equations will be expressed in a dimensionless form, by introducing 
the following scaled time, velocity, and displacement variables: 

t 

T = -, 

T 

04, vg = —{v x ,v y ), (ii) 

gr 

(X,Y) = -^(x,y), 
so that Eq. (flOl) leads to the dimensionless system 

V x = -v x , 



(12) 



Vy = -1-V y , 

which, once inserted into the dimensionless counterpart of Eq. 

d 2 Y 1 d (Vy\ _ VyV X ~ V X Vy 



dX 2 V x dT \ V X J Vi 

leads to the following differential equation for the trajectory: 

d 2 Y (-l~V y )V x + V x V y J_ 
dX 2 K 3 V 2 ' 



(13) 



(14) 



formally identical to Eq. ([7]). However, differently from the previous case, the 
x-component of the velocity is no longer constant, and we have to find the 
dependance of V x on the variable X. This can be easily done on using the 
first row of Eq. (fl"2|) together with the chain rule in Eq. ([5]) to remove T, thus 
obtaining 

t - ™ 

which gives 

V x = C - X , (16) 



3 This is equivalent to formally set r = 1 and g = 1, in Eq. {Tojl . 
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where the constant C can be found by imposing the initial condition on V x , i.e., 
V r x (0) = Fo, x , so that 

V X (X) = V , x - X. (17) 
On substituting from Eq. (fTTj) into Eq. (|14p we obtain 

dX 2 (Vb ia - x) 2 ' 1 0J 

that can be solved by using elementary quadratures. In particular, on taking 
the initial conditions in Eq. ([8]) into account, after some algebra we have 

Y = Xt!m6 + +logfl--^-V (19) 

which, on using Eq. (|llj) in order to remove the normalization and to come back 
to physical quantities, becomes 



y(x) = x tan# + gr 2 



loe 



(20) 



_V0,xT \ V 0}X 

It should be noted that Eq. (|20p is known from the literature, although its 
derivation is customarily carried out by first solving the time-domain dynamical 
equation (|10[) and only after removing the variable t among the functions x = 
x(t) and y — y(f)0 It is also worth showing to the student how, for small values 
of drag force, the shape of the trajectory tends to become identical to the ideal 
free- friction case of Eq. ©. To this end it is sufficient to consider the limiting 
expression, for r — > oo, of the logarithmic function in Eq. (|20[) which, by using 
a second-order Taylor expansion, turns out to be 

£ + log(l-0 ~ ^0. (21) 



4 Trajectory with a quadratic-speed drag force 
4.1 Derivation of the differential equation 

When the resistant medium exerts on the body a force proportional to the 
squared modulus of the speed, Newton's law takes on the following form0 

a = g - - vv , (22) 



4 See for instance Refs. 181123). 

5 It could be worth putting into evidence the fact that, as reported for instance in Ref. 1241 . 
Newton himself was aware that a linear resistance model was not the most suitable from a 
physical viewpoint, but rather the resistance offered by media without rigidity would have 
been expected to be proportional to the square of speed, as he wrote at the end of Sec. I of 
Book II of his Principia |25j : 

Cceterum resistentiam corporum esse in ratione velocitatis, Hypothesis est 
magis Mathematica quam Naturalis. Obtinet hcec ratio quamproxime ubi cor- 
pora in Mediis rigore aliquo pr&ditis tardissime moventur. In Mediis autera quee 
rigore omni vacant resistentieE corporum sunt in duplicata ratione velocitatum. 



4 Trajectory with a quadratic-speed drag force 



7 



where I is a parameter whose physical dimensions are those of a length. Similarly 
as done for the linear resistant medium, it is worth giving Newton's law a 
dimensionless dress by introducing the following scaled variables: 

(V x ,Vy) = -±=(V x ,Vy), (23) 

V5« 



(X,Y) = \(x,y), 
so that Eq. (|2"2"j) takes on the form [IT] 

V x = -VV X , 



(24) 



Vy = -1~ VVy , 



where the speed V, defined as 



V = JVi + V* , (25) 



acts as a coupling factor. On substituting from Eq. (|24l) into Eq. (ITS")) we find 
again 

d 2 Y _ {-l-VV y )V x + VV X V V __ 1 



dX 2 V$ V 2 ' 



(26) 



but now it is no longer possible to achieve an explicit X-dependance for V x , a 
fact which prevents Eq. (|26l) to be solvable via simple quadratures. In fact, on 
applying the chain rule to the first row of Eq. (f24|) we have 

dV T dV r 



,i7 Vx di = ' VxV = ~ Vx ^ + v * ' (27) 



so that 



dV x 



(1V -jVi + V*. (28) 

Moreover, on deriving both sides of Eq. (j2l)]) with respect to X, on taking 
Eq. (|2"5|) into account, so that 



d3r = 4^ = -4 Ji + r^v = -4 Ji + r^v, (^) 



dx 3 v; 3 dx v 2 V V 

and on using again Eq. (|26p we obtain 



^ 2S,l + « 2 , (30) 



dx 3 dx 2 V V d ^ / ' 
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which is the differential equation satisfied by the Cartesian representation of 
the body trajectory Y = Y(X). It has to be solved together with the initial 
conditions 

r(o) = o, 

Y'(O)=tan0 o , (n) 

v 0,x 

but unfortunately the solution of the Cauchy problem in Eqs. ([30]) and (131|) 
cannot be given a closed form; nevertheless, nowadays it is easily achievable 
numerically, due to the availability of several commercial computational plat- 
forms. 

From a didactical point of view it could be worth showing how Eq. (|30|) leads 
in elementary way to a well known analytical approximation of Y(X), called by 
Parker the "short-time approximation" [TTJ @ which is valid in the limit of small 
velocity slopes. In fact, within such approximation it is sufficient to neglect the 
squared term into the square root of Eq. (1301) , in such a way that the differential 
equation reduces to 

^~2^ (32) 
dX 3 ~ dX 2 ' 

which can be solved simply by iterated elementary quadratures. In particular, 
on taking Eq. (1311) into account, after simple algebra we obtain the approximated 
solution, say Y s t, as 

which coincides with the equation given at the end of Sec. Ill in Ref. [TTJ and 
with Eq. (26) at p. 297 of Ref. [23]. 

4.2 Power series expansion of the trajectory 

It must be appreciated that the mathematical structure of Eq. (I30[) turns out 
to be particularly suitable to derive in a simple way the Taylor series expansion 
of the trajectory, namely 

oo 

n*) = £f^ fe , (34) 

fc=0 

where = Y^ (0) denotes the fcth-order spatial derivative of Y(X), evaluated 
at X = 0. In fact, starting from the initial conditions given in Eq. (I3TT) we have 



6 As the name suggests, the derivation of Parker was done in the time domain. A different 
derivation can also be found in the beautiful book by Lamb |23l . 
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for the coefficient the expression 



a 3 = Y"'(0) = 2Y"(0) v^TFW = -77- 2 T , (35) 

v 0,x cos fo 

which differs from the corresponding coefficient in the power series expansion of 
Ygt by a factor cos0o that, in the case of small slopes, can be replaced by the 
unity. As far as the coefficient is concerned, it is sufficient to derive both 
sides of Eq. (f30|) with respect to X to obtain 

VI + Y> 2 

which, after rearranging and simplifying, gives at once 

A f o sin 0n \ 

whereas the short-time approximation would provide 

^(0)=-4-- (38) 

v 0.x 

It is clear that, in principle, all terms of the power series expansion in Eq. (|34l) 
could be evaluated by iterating this approach; however, on increasing the trun- 
cation of the expansion the complexity of the analytical expressions of the co- 
efficients afc rapidly grows0 although it is easy to implement their evaluation, 
up to arbitrarily high orders, by using any commercial symbolic computational 
platform. 

To give an idea of the performances of the above polynomial approximation, 
we analyze a case already studied in the past, related to the motion of a volley- 
ball in air [5 . Figure [2] shows the ball trajectory (within dimensionless spatial 
variables) for an initial velocity with components (vo tX ,Vo tV ) = (2,5) m/s in a 
medium characterized by £ ~ 27.64 m. Such data refer to the results shown in 
Fig. 4 of Ref [5]. The exact solution, obtained by solving numerically Eq. (|30[) 
through the standard command NDSolve of the symbolic language Mathematica 
is represented by the circles, while the dashed and the solid curves are the poly- 
nomial approximations obtained by truncating the series in Eq. (|34j) up to k = 3 
and k = 4, respectively. The dotted curve represents the trajectory correspond- 
ing to the ideal case in vacuum. From the figure it is seen that the fourth-order 
polynomial displays an excellent agreement with the exact curve along the whole 
range of interest of a;. To make the comparison between our figure and Fig. 4 



7 As a further example we give the expression of the 5th-order coefficient: 
8 cos6>o „ 1/2 ^ a , f 



a 5 = I Vq - 2V<f sin O + - cos" 8 



which should be compared to its "short-time" version, given by — 8/Vg 2 , 
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Fig. 2: Trajectory (within dimensionless spatial variables) of a volleyball in air 
(£ ~ 27.64 m) as numerically evaluated from Eq. (|30|) (open circles) to- 
gether with the 3rd-order (dashed curve) and the 4th-order (solid curve) 
polynomial approximation obtained by truncation of the power series 
expansion in Eq. (|34j) . The upper curve refer to the initial conditions 
(vo, xi v o,y) — (2 ; 5) m/s, while the other curve to the initial conditions 
( v o,x, Wq, j/) = (6, 3/2) m/s, in order to be comparable to Fig. 4 of Ref. [5]. 
The dotted curve represents the trajectory corresponding to the ideal 
case in vacuum. 

of Ref. 5 more complete, the trajectory corresponding to the initial condition 
( v o,x, v o,y) = (6, 3/2) m/s is also plotted. In this case the 3rd- and the 4th-order 
approximants are practically undistinguishable. 

Figure [3] shows the results corresponding to up = 14 m/s and to several 
values of the initial angle 8q, retrieved from Fig. 6 of Ref. 5 . In particular, 
for each curve it is shown the exact solution (open circles) together with an 
iVth-order polynomial approximation (solid curve), with N (labelled for each 
curve) being the lowest value of the truncation order of the power series which 
guarantees a satisfactory agreement (at least at a visual level) to be achieved. 

4.3 An elementary derivation of the hodograph 

In the present section we want to show how the third-order differential equa- 
tion (|3T?|) leads in a simple way to the hodograph of the motion [THl HO] • To this 
end we first recast Eq. (|3"0)) in terms of the slope i] = Y' as follows: 



(39) 
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Fig. 3: The same as in Fig. [2j but for vo = 14 m/s and for several values of 
the initial angle 9o, extracted from Fig. 6 of Ref. [5]. for each curve 
it is shown the exact solution (open circles) together with an iVth-order 
polynomial approximation (solid curve) , with N (labelled for each curve) 
being the truncation order of the power series which has been chosen in 
order to allows a satisfactory agreement to be achieved. 



which, on letting 

(40 

becomes 

dry' 



iY 277' Vl + V 2 - (41) 

Although the differential equation cannot be solved in analytical terms, it is still 
possible to find an explicit functional expression for the quantity 77' whether it 
is thought of as a function of 77. Again the trick is to use the chain rule, by 
letting 

_d_ = d77_ _d_ = _d_ 

dX dX dry ' dry ' V ' 



dT?' 



so that Eq. (|4Tj) becomes 



-^ = 2v/rT^, (43) 
Q77 
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and rf can be retrieved through a simple quadrature, which gives 
r\ = rj y/l + rj 2 + arcsinhfi + C = 



(44) 



= j] y/1 + rf + log(r? + y/i + rf) + C . 



The constant C has to be determined by imposing that at the start of the 
trajectory, when its slope is 770 = tan0n, the corresponding value of 77', say r/' 
is, according to Eq. $Zjo§, 

^ = -4 = ~*^v (45) 

so that, after simple algebra, we obtain 

C = -^r-f(9 ), (46) 

v Q.x 

where the function /(•) is defined by 

„,„, sin0 , l + sin# , 

= —2H + lo S IT ■ (47) 

cos cos 

On substituting from Eqs. (f46|) and (|47l) into Eq. (|44|) the following differential 
equation for the function r/(X) is then obtained: 

|| = + V y/l + V 2 + log(?7 + V^Tv 1 ) - f(9 ) , (48) 

which, on recalling that 77 = tan 9 and that the l.h.s. does coincide with Y"(X), 
after algebra can be recast as follows: 

S=-4 +/w - /( *° } ' (49) 

and, on taking Eq. (|26p into account, after rearranging leads to the well known 
expression of the hodograph [26 

] — — ■ (50) 



cos9^l + V 2 x [f(9 ) - f(9)} 



5 Conclusions 



An elementary exposition of the trajectory determination problem for a body 
launched at an arbitrary direction and subject to the gravity in the presence 
of a resistant medium has been proposed. The main task of such exposition 
was to limit at most the mathematical complexity in order to be appreciable 
also by first-year undergraduate students, especially in the case of a quadratic 
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resistant medium. To this end, the use of time domain Newton's laws of the 
motion was systematically avoided by removing, through the use of the deriva- 
tive 'chain rule', the temporal variable before solving them. In this way the 
differential equation for the representation of the trajectory has been derived in 
a Cartesian reference frame which turns out to be naturally suited to represent 
the trajectory shape. In the study of the motion in vacuum and in a resistant 
medium linear in the velocity the proposed approach was able to retrieve, in an 
elementary way, the well known exact expressions of the trajectory equation. 
For a body moving in a resistant medium which exerts a force quadratic with 
respect to the velocity, the proposed approach allowed the differential equation 
for the Cartesian representation of the trajectory to be derived in a rather sim- 
ple way. From the same differential equation simple polynomial approximants 
of the trajectory, having the form of truncated power series expansions, have 
been built up to arbitrary values of the truncation order and compared to the 
results recently obtained through a homotopy analysis based approach. Finally, 
to complete such a sort of changeover toward more advanced mathematical 
treatments of the problem, our approach was also used to derive, again in an 
elementary way, the closed form of the motion hodograph, a well known result 
from the literature. 
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